Tack energy and switchable adhesion of liquid crystal elastomers 
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The mechanical properties of liquid crystal elastomers (LCEs) make them suitable candidates 
for pressure-sensitive adhesives (PSAs). Using the nematic dumbbell constitutive model, and the 
block model of PSAs, we study their tack energy and the debonding process as could be measured 
experimentally in the probe-tack test. To investigate their performance as switchable PSAs we 
compare the tack energy for the director aligned parallel, and perpendicular to the substrate normal, 
and for the isotropic state. We find that the tack energy is larger in the parallel alignment than the 
isotropic case by over a factor of two. The tack energy for the perpendicular alignment can be 50% 
less than the isotropic case. We propose a mechanism for reversibly switchable adhesion based on 
the reversibility of the isotropic to nematic transition. Finally we consider the influence of several 
material parameters that could be used to tune the stress-strain response. 



I. INTRODUCTION 

Pressure sensitive adhesives (PSAs) are soft materials 
that adhere to nearly any surface when low pressure is 
applied. Their effectiveness can be described by the total 
energy required to separate the adhesive from a surface, 
known as the tack energy. Experimentally this can be 
measured in the probe-tack test, where the force required 
to remove a probe moving at constant velocity from 
an adhesive film is measured and a stress-deformation 
curve is produced. An experimentally measured stress- 
deformation curve produced from a probe-tack test on 
two materials is shown in fig. [T] For small deformations 
the force rises rapidly with extension up to a peak force. 
During this phase cavities form within the adhesive and 
grow. Following the peak force there is a pronounced 
plateau in the force which can be upward pointing (PI) 
if the material strain-hardens. During this phase fibrilla- 
tion occurs and the force is dominated by the viscoelas- 
tic properties of the adhesive rather than the cavities. 
Ultimately the plateau region ends when the material 
detaches from the probe. Detachment can occur at the 
surfaces, which is known as adhesive failure and occurs 
for material PI in fig. [TJ Alternatively detachment can 
happen in the bulk, which is known as cohesive failure 
and occurs for material P2 in fig. [T] 

To achieve a high tack energy the PSA must have a 
low dynamic modulus (typically 0.1 MPa at 1 Hz) to 
make conformal contact with the substrate @|. As the 
probe is retracted, the adhesive film is drawn into fib- 
rils. The adhesive must be sufficiently soft so the fibrils 
do not detach at small strain, yet still require some en- 
ergy to deform Q. PSAs are typically made from high 
molecular weight polymers that are lightly crosslinked to 
form viscoclastic solids. Acrylic, styrcnic, and siloxane 
based polymers have been refined for use as PSAs [||. 
The bulk rheological properties of a PSA are important 
in determining its tack energy. The optimal stress-strain 
behaviour is initially strain softening, to aid crack blunt- 
ing, and eventually becomes strain hardening to stiffen 
the fibrils, and increase the force required in the latter 
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FIG. 1. The stress-deformation curve measured by a probe- 
tack test performed on two materials, PI: Poly(butyl acry- 
late) copolymer and P2: Poly(butyl acrylate) copolymer with 
CTA. Probe-tack measurements were performed at constant 
retraction velocity of 0.1mm s _1 [lj. 



stages of debonding 

Gay et al. have presented a successful theoretical de- 
scription of PSA debonding based on the growth of bub- 
bles at the substrate/adhesive interface [fj]. These bub- 
bles join, forming larger bubbles which are then drawn 
into fibrils. The main physical processes of homogeneous 
deformation, nucleation of cavities, followed by fibril for- 
mation have been included in a simplified model of the 
debonding process called the Block model @, [!]■ Here 
the adhesive layer is divided up into N equal rectangu- 
lar blocks that can undergo a combination of shear and 
stretching deformation, as well as slipping on the sub- 
strate. The force-displacement curves of the block model 
reproduce the characteristic features seen in probe-tack 
experiments. 

There is considerable interest in being able to turn off 
the adhesive properties of PSAs using external stimuli 
such as light, humidity, and temperature. For example 
methacrylate-functionalizcd adhesives containing a pho- 
toinitiator can show an almost complete loss of adhe- 
sion when irradiated under a halogen lamp, as a result 
of photo-initiated crosslinking raising the elastic mod- 
ulus 0, [poj . The switching here however is only one- 
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way, the adhesive cannot be returned to a tacky state. 
Trenor et. al [TTJ achieved two-way switching by us- 
ing courmarin-functionalized acrylatc adhesives, whereby 
UV-A radiation was used to switch off adhesion and UV- 
C radiation was used to (partially) switch it back on. 
Altering ambient humidity has been used to reversibly 
modify the surface composition of polymer blends thus 
changing the adhesive strength [l2|; note that here the 
adhesion is modified by altering the substrate rather than 
the adhesive itself. There are several other examples 
of switchable adhesion/ wetting based on polymer brush 
surfaces driven by an external stimulus [13l . Il4| . Liquid 
Crystal Polymer (LCP) based adhesives which undergo a 
reversible Smectic-Isotropic transition have been shown 
to have a transition between a tacky and a non-tacky 
regime as one cools to the smectic phase [l5j . The low 
temperature smectic phase is harder, and less wet ting 
than the isotropic phase. Kamperman and Synytska [16| 
provide a review of two mechanisms to achieve switchable 
adhesion; topography and chemical functionality. 

In this paper we propose an application for weakly 
crosslinked LCPs, called liquid crystal elastomers 
(LCEs), as anisotropic adhesives whose tack energy can 
be switched depending on the orientation of the director, 
and the degree of liquid crystalline order. The mecha- 
nism here is based on the change in the bulk rhcology of 
the adhesive. 

LCEs are unique materials that couple liquid crystal 
mesogens to the underlying polymer network. In the high 
temperature isotropic phase, they behave much like con- 
ventional rubbers. When cooled into the nematic phase 
the orientation of the director is crucial to determining 
their mechanical behaviour. When stretched parallel to 
the director they behave like uniaxial solids. On stretch- 
ing perpendicular to the director they exhibit (semi-)soft 
elasticity - that is strain increases with little increase 
in stress. They have a long plateau in their stress-strain 
curve before undergoing strain hardening to values of the 
modulus found in isotropic rubbers at low extension. The 
nematic_phase has been synthesized using both acrylate 
based [lTj and siloxane based polymer chemistry fl8j . 
though in practice the siloxane based LCEs have a lower 
T g and more durable mechanical properties. Their me- 
chanical behaviour is well described by a phantom chain 
model of a network of anisotropic Gaussian chains fl9j . 
This description predicts the plateau in the stress-strain 
curve associated with soft elasticity 20]. However this 
equilibrium statistical physics model does not describe 
the dynamics of the LCE that are essential to model ad- 
hesive properties. We will show that the nematic dumb- 
bell model of Maffettone and Marrucci includes the equi- 
librium behaviour of nematic elastomers, and also pro- 
vides a description of their dynamics [2l| . By using this 
constitutive model in the block model we describe the 
behaviour of LCE based viscoclastic adhesives. 



II. MODEL SYSTEM 

Here we consider an idealised version of the probe-tack 
test. A thin adhesive layer is placed between two rigid 
substrates, which are then pulled apart with a constant 
velocity. The force required to separate the substrates 
is measured as a function of the displacement between 
them. We will model the subsequent debonding process 
in two dimensions by assuming that the adhesive and any 
deformational flow are confined to the xy plane. The 
debonding process involves several complex phenomena, 
but the following dominant effects will be included in 
our description; large strain deformation of the adhesive, 
cavity expansion/contraction, slip of the adhesives at the 
substrate interface. These effects have been incorporated 
in the Block model developed by Yamaguchi et al. 0], 
which gives a 2-dimensional description of isotropic ad- 
hesives at low Reynolds number. We modify the Block 
model by utilising simpler shapes and deformations for 
the blocks, and a material constitutive law that describes 
LCEs and their associated director reorientation. 



A. Block Model 

The initial adhesive layer has height Hq and length 
Lq (we assume the adhesive is thin, i.e. Lq » Hq). 
This layer is divided into TV rectangular blocks of equal 
width Wq = Lq/N, which can undergo slip at the inter- 
face with the substrate. Each block is assumed to un- 
dergo a stretching deformation along the y-direction, no 
deformation in the z-direction, a corresponding volume- 
conserving contraction in the x-direction, and a piece- 
wise linear simple shear deformation in the x-direction as 
shown in Fig. [2j The motion of the ith block is charac- 
terised by three parameters: the elongation along the y- 
direction denoted A, the shear deformation denoted \ xy ,i 
and the x-coordinate of the centre of mass of the block 
denoted Xi. 

On stretching the adhesive layer by a factor A in the 
y direction, each block has a new height H = XHq and 
a new width W = Wq/X. The shear deformation of a 
block can be described by dividing it into two sub-blocks 
of width W and height H/2. The lower sub-block shears 
by X X y while the upper sub-block has the opposite shear 
— X xy as illustrated in Fig. [5] This deformation is slightly 
different to that employed by Yamaguchi et al. [3] • They 
allowed a parabolic deformation in the x-direction, which 
implies a shear which varies continuously with position 
along the y-axis. We anticipate a coupling between shears 
and the director in our nematic adhesives, and modelling 
the nematic degrees of freedom is substantially simplified 
if we assume that the sub-blocks have constant shear and 
director. We will see later that this change simply mod- 
ifies the constant prefactor in the equation for the shear 
of the block. The deformation gradient tensor A for each 
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FIG. 2. The deformation of an individual block involves an 
area preserving elongation along the y direction, and a piece- 
wise linear deformation in the xy plane. The x component of 
the centre of mass of the block is Xi while the point of contact 
at the interface between the block and the substrate is X B i. 
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FIG. 3. (a) If the difference (Xi+i — Xi) is greater than Wo/A 
we say a cavity has formed between the blocks as illustrated. 
The pressure Pi on the faces of blocks i and i + 1 is required 
to calculate the debonding force, (b) We use a simplified 
model of the cavity, assuming it to be a spherical bubble of 
radius Ri placed in an infinite Newtonian fluid with viscosity 
r\ — Gt. Assuming creep flow and applying force balance at 
the interface between the bubble and the interface leads to 
a value for Pi, the pressure on the fluid side of the bubble 
interface. 



sub-block can be written as 

I ±A : 



A = 



A 



(1) 



where the sign of the shear X xy depends upon which sub- 
block is being described. 

Let (X S i,0) be the mid-point of the contact between 
the block and the substrate in the deformed state, and 
(x s i , 0) the equivalent point in the undeformed state. 
Then transforming between the undeformed state (x, y) 
and deformed state (X,Y) using (dX,dY) = A • (dx,dy) 
and integrating we obtain 

X sl + + X xy . iy y < H /2 

^ X sl + + X xy ,i(H -y) y> H /2 

Y(x,y) = Xy. 

Hence the cc-coordinate of the centre of mass of the ith 
block, is given by 

ry=H i-x=x si + W /2 

Xi= / X{x,y)dxdy/(H Wo) 

Jy=0 Jx=x si -Wo/2 



Xi(x,y) 



v I \ ^0 



(2) 



Finally, we will require the elements of the velocity gradi- 
ent tensor in the deformed configuration Kij = \7jVi. In 
terms of the deformation gradient tensor A this is given 
by K = A • A -1 : 
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B. Cavity Expansion 



(3) 



As the separation between substrates XH increases, 
there is a corresponding reduction in block width Wo/ A. 



If the displacement between the centres of the blocks 
(Xi+i — Xi) is not equal to their width, a gap is cre- 
ated between the blocks labelled i + 1 and i. We treat 
this gap as a circular cavity with radius Ri having the 
same area as the void between the blocks, see Fig. [3] The 
instantaneous radius of the circular bubble is given by 



H a X I X 



X,. 



Wo 
A 



(4) 



The Rayleigh-Plesset (RP) equation [22|, [23| describes 
the dynamics of a spherical bubble placed in a Newtonian 
fluid with viscosity r\. If Pi is the pressure in the fluid 
at the bubble interface, P cav ,i the pressure in the cavity 
at the fluid interface and 7 the surface tension, then the 
RP equation is 



dRi _ -R* / _ _ J_ 

at 2r) \ Ri 



(5) 



If the cavities do not contain air then P cav ,i is simply zero. 
If the cavity includes air, initially at atmospheric pressure 
P atm in the undeformed state then P ca v,i is related to the 
area change in the cavity via 



Pcav,i(t) — -^a 



-Ri(O) 
Ri(t) 



(6) 



This is a highly simplified model of the 
growth/contraction of the cavities, it ignores both 
the elasticity and the anisotropy of the viscoclastic 
adhesive surrounding the cavity. Given the anisotropy 
of the adhesive it is unlikely bubbles would grow in a 
uniform circular/spherical way. Nevertheless the model 
has the correct qualitative features for any model of 
cavity dynamics and we adopt it for simplicity. 
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C. Block positions 

Equation (|4|) gives the radii of the N — 1 cavities in 
terms of the positions of N blocks. To invert this equa- 
tion and hence determine the positions of the N blocks 
we need an additional constraint. The applied external 
stretching force is in the y direction, there is no external 
force in the x direction. The x coordinate of the centre 
of mass of the N blocks thus does not change with time 



N 



N 



'£x i (t) = '£x i (t = o) = o. 



(7) 



i=l 



i=l 



Therefore if we know the cavity radii, we can use equa- 
tions (@| and([7]) to solve for the positions of all N blocks. 



D. Stress Tensor 

The total stress tensor within each block is made up of 
two components, the polymer stress arising from forces 
transmitted by the polymer chains within the adhesive 
denoted by £, and the isotropic pressure term —pi. The 
total stress tensor <r is produced by adding these two 
components 



(7 = - P I + S. 



(8) 



The constitutive equation obeyed by S will be discussed 
in §1111 We assume that the components of the polymer 
stress tensor are homogenous within each block. This is 
not so for the pressure contribution. The ith block has 
two free surfaces. At these surfaces the xx component of 
the total stress tensor is a xx = —Pi and a xx = — P»+i. 
To accommodate this change in the total stress it is clear 
the term —pi must vary across the block. We will assume 
this variation is linear. 



E. Slippage at the interface 

We denote by a S i the shear stress at the interface be- 
tween the ith block and the substrate. We assume a 
simple linear relationship between the position of the in- 
terface between the block and the substrate X S i and the 
shear stress <r si at the same place 



d Y 



(9) 



Recall X s i = Xi — X xy ^Ho/A, thus having previously de- 
termined the positions of the blocks Xi we can rewrite 
this as an equation for the rate of change of the shear 
strain A^, i . 



F. Force Balance 

The difference in the pressure on cither side of the 
block (Pi+i — Pi) is balanced by the shear stress, leading 



to a force balance equation 

(P i+1 - P,)H Q X 



-2a. 



Wo 
; A ■ 



(10) 



Once the shear stress of each block has been calculated, 
we can use this equation, coupled with the boundary con- 
dition Pi = P a tm to calculate the pressure within each 
block. 



G. Debonding Force 

To determine the debonding force F y (the total force 
acting on the substrate) we require the pressure, and the 
xx and yy components of the polymer stress tensor S. 
The a xx component at the left and right edges of the ith 
block can be used to calculate an approximate value of 
the average pressure in the block p using the assumption 
of a linear variation in p through the block 



\x si -W /2 



P 



x Bl -W /2 



" xx > l \x si + W /2 ~ " xx > z P \x 3 ,+W a /2 

P = P\x 3i -W /2 + P \x si + W /2 

P + Pj-1 
~ ^xx,i t 2 

The (Jyy.i component from the ith block is thus 

P. + Pi-, 



Pi-i (11) 
P (12) 



(13) 



(14) 



This result contains the normal stress difference of S, 
hence the addition of an isotropic pressure term to the 
constitutive model will be absorbed into the pressure p 
of Eq. (fT3")) . The total force F y can be calculated by 
summing the force due to atmospheric pressure on the 
substrates, and that due to each block in the adhesive 



A 

Fy = ~Y 



Pi + Pi 



i+i 



i=l ^ 

(15) 

where A is the initial contact area. Note Yamaguchi et 
al. did not use a normal stress difference in Q, however 
it makes little difference to the results for the parameters 
and strain ranges used there. Dividing Eq. (fTS"]) through 
by A we obtain the engineering stress in debonding 



f„„ = — ^- = i^yy ^Jxx) ~\~ (-Patm 



yy 



P), 



(16) 



where T represents an average of quantity T over all 
blocks. To calculate the tack energy £t we integrate the 
debonding force over the distance moved by the upper 
substrate 

£t = f V Fydy' = AH Q f a(X')d\' . (17) 

Jy' = H a J\' = l 
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III. CONSTITUTIVE RELATIONSHIP 

To determine the tack energy of the PSA we require the 
rheological properties of the adhesive layer in mi Al As we 
are modelling a LCE based PSA, we will use the nematic 
dumbbell model of Maffettonc and Marrucci [2l|. We 
consider a polymer chain made up of TV freely hinged ne- 
matogcnic rods with individual length b. The end-to-end 
vector for the chain is R. Assuming Gaussian fluctua- 
tions the quantity of interest is the scaled second moment 
of the end-to-end vector, that is W = 3 (RR) /Nb 2 . The 
dynamics of W arc governed by 

w = -i -!- (r 1 ■ w + w-r 1 ) , (is) 

H = G£~ 1 -W, (19) 

v 

where W= dW/dt — K • W — W • K is the upper con- 
verted Maxwell derivative (UCMD) and K is the velocity 
gradient tensor in Eq. ([3]) . The inverse chain shape tensor 
(T 1 is given by 

< - 1 -<rbj( I -iHs fifi )' <20) 

in which the nematic director n describes the average 
orientation of the nematogenic units and the order pa- 
rameter S which describes the degree of alignment along 
the director (S = 1 corresponds to perfect ordering, while 
5 = corresponds to an isotropic phase). The polymer 
stress X is that arising from the polymer chains within 
the material. The time-scale r is related to the diffusivity 
of the chain ends D via r = Nb 2 /6D, while the modulus 

G = ck B T, (21) 

where c is the number of chains per unit volume. 

This Gaussian model does not include effects of finite 
chain extensibility, or entanglements. However Gaussian 
models have been shown to provide a good description of 
the mechanical behaviour of liquid crystalline elastomers 



Any isotropic stress can be added to the definition of 
S, as it can be absorbed into the pressure term in 
Eq. (jHJ. It is convenient to subtract GI from the stress 
£' = £ — GI. This subtraction is useful since £' is zero 
in the undeformed state. 

v , £' 

£' +— = 2GD, (25) 

T 

where D = (K + K T )/2 is the symmetric part of the 
velocity gradient tensor. This is identical to the consti- 
tutive relationship used by Yamaguchi et al. Q . 

B. Director dynamics 

To complete our description of the dynamics of the 
polymer stress we must also describe the behaviour of 
the nematic director with time. In general this is a com- 
plicated task involving the coupling of the director to 
the flow field. Maffetone and Marrucci [2l[ identify two 
simplified regimes 

1. Weak external field. In this case the polymer stress 
tensor is required to be symmetric. To achieve this 
in Eq. (|T9"|) the director n must be one of the eigen- 
vectors of W. Since in static equilibrium we re- 
quire the polymer stress to be isotropic, W eg must 
be equal to £, from which we can identify n as the 
eigenvector corresponding to the largest eigenvalue 
of W. 

2. Strong external field. In this case the torques aris- 
ing from coupling to the flow field are insufficient to 
move n from the direction imposed by the external 
field. The polymer stress will not be symmetric in 
this case. 

We model the first of these limits, where the director 
responds much faster than the polymers. We leave the 
more complicated task of generalised director dynamics 
to future work. 



A. Isotropic Limit 

Yamaguchi et al. first described the block model using 
the constitutive equation of an isotropic Maxwell fluid for 
the polymer stress. Here we demonstrate the isotropic 
limit of Eqs. (fT5| and (fT9"f produce the same constitutive 
equation. 

To obtain the isotropic limit, we set S = in Eq. (fT8|) 
and ([15]). This gives 



v 

W 



V-W, 

T T 

GW. 



(22) 

£ = GW. (23) 
Substituting for W in equation (f22j) in terms of £ 



v 
£= 



G, 



1. 



(24) 



C. Quasi-Static Limit 

There are consistent theoretical descriptions of nematic 
elastomers derived phcnomcnologically from continuum 
mechanics [24|, and from a microscopic equilibrium sta- 
tistical physics model [lj|. These models both contain 
the Goldstone modes predicted in Nematic elastomers us- 
ing symmetry arguments [25| , known as soft modes. The 
microscopic model produces the following trace formula 
for the free energy density F 



F = ±GTr 



X- £ ■ X t 



(26) 



where G is the shear modulus of the rubber defined in 
Eq. (|2ip , £q is the initial chain shape distribution, and it 
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is assumed that det[A] = 1. The soft modes permitted 
by this free energy have the following explicit form 



A = £2 . Q .£ Q 



(27) 



where Q is an orthogonal tensor. They arise because 
states which are related to each other by a simple rotation 
of the chain shape distribution tensor £ have the same 
energy. The true stress can be derived from Eq. ([26|l by 
differentiating F with respect to A, then post-multiplying 
by A T , producing a polymer stress component 



GC 1 ■ A -£o • A T . 



(28) 



Eqs. (|18jl and (|T9|) also permit soft mode solutions. Mo- 
tivated by Eq. (|2"5|) we substitute W = A • £ ■ A T , into 
Eq. (jT5)) and observe that the UCMD is identically zero, 
leaving 







A-^ -A J +X-£ -X 1 ■£-') 
2t 



(29) 



This equation holds provided that A obeys Eq. (|27|). 
The stress tensor associated with these modes is simply 
S = GI. The normal stress difference Yj vv — E xx = 0, 
and thus these modes would not contribute directly to the 
debonding force of Eq. (JTSJ) (they contribute indirectly in 
so much as the pressure difference across the block is de- 
termined by the shear stress associated with the mode). 
In typical experiments on nematic elastomers stretching 
perpendicular to the nematic director produces a soft re- 
sponse as the director rotates towards the stretch direc- 
tion. Stretching parallel to the director produces a hard 
elastic response. We thus have reason to suspect that 
the debonding force could be substantially different for 
parallel and perpendicular geometries. 



IV. NUMERICAL METHOD 

To solve our set of differential equations we implement 
a mixed implicit /explicit finite difference scheme to step 
forwards in time. To improve computational accuracy 
we reduce the number of free parameters by combining 
those with equivalent effects. The numerical results are 
in agreement with the analytic single block results and 
the semi-analytic two block case with t — > 00. They are 
consistent with the isotropic results of Yamaguchi et al. 



Initial Condition 



The initial condition has N blocks all of which have 
A = 1 and X xy .i = 0. There are N — 1 cavities with ini- 
tial radius Ri = Rq exp(xi) where Xi ' 1S assigned from a 
normal distribution with (%) = and (x 2 ) = 1 an( i ^0 
is a characteristic size. The initial positions X,(i = 0) 
can then be determined by solving Eqs. ^ and (0. 



The pressures at the block interfaces are initially set 
to atmospheric pressure, thus Pj = P atm - We assume 
that the cavities form at the substrate/adhesive inter- 
face and are initially filled with air, and thus we also 
set the cavity pressure equal to atmospheric pressure, 
Pcav,i = P&tm [26| • Experimentally, the cavities are some- 
times found to be filled with vacuum even if the cavities 
formed at the substrate/adhesive interface [27j. Yam- 
aguchi et. al found the presence or absence of air within 
the cavities made very little difference to the resulting 
tack curves Q. We divide the time t into units St and 
advance the equation set using a mixed explicit/implicit 
scheme. We introduce a superscript to quantities to in- 
dicate at which time-step they are evaluated, i.e. A^™' is 
evaluated at the nth time step when t = nSt. 



1. We increment 

A (n+1) = A (n) , 



Time Stepping 

the strain 



explicitly 



2. We update the cavity radius semi-implicitly 



R 



(n+l) 



2r)R 



>(") 



3. Using the new values of the cavity radii we obtain 
the cavity pressure 



>0) 



R 



(») 



R 



(n+l) 



4. With the new values for the strain A(" +1 ' and cavity 
radii {R^ +1 } we solve Eqs. (2]) and ([7]) for the new 
positions of the blocks {X^ n+1 ^}. 

5. Using X si = Xi - A, 
Eq. © 



A 



(n+l) 
*xy,i 



A 



(n) 



4 



xy,%Ho/4: we update A^j from 

(n) 



(n+l) 



a: 



1 1 



-St 



This equation replaces the curvature of the blocks 
used by Yamaguchi et al.. It describes the same 
physical process of block shear, though it differs by 
a geometrical factor. 

6. With the shear strains {Xxy +1 ^} and A( n+1 ) we cal- 
culate the elements of the deformation gradient ten- 
sor 



K 
K 

jX 



A 



A(™+!) 
-K, 



yy 
K, 



yy 



(A 



(n+l) 
xy,i 



xy,i J 



\( n +V5t 
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7. We now update elements of W using a mixed ex- 
plicit/implicit scheme 



Parameter (symbol) 



Value 



yy 



W(«+i) 

' 1 xy 



ir, 



(«) , St 



St „( n )A n )vtr( n ) 



l-2K yy St+^ 



i i (™) (™) 
1 + an y n y 



1 + 



Si 



2(1-S) 



1 + (i + 1) n x n) n x n) 
n x n) n y n) (W^ +W^ +1) ) 



2(1-5) 



+5tK xy W^ 



W/0+ 1 ) 



1 + 2KyySt + 



, (n) (n) 

1 + cmi 'ni 



11 3 1-5 



+2<^ y W^ +1 ) 



where a 



3S 
2S+1 ■ 



8. Determine the director by finding the eigenvector 
of \y(" +1 ) corresponding to the largest eigenvalue. 

9. Calculate the updated polymer stress from 
Eq. PI). 

10. Calculate the new pressure field {P- n+1) } from 
Eq. (p|) with the boundary condition Pi = P a t m - 

11. Repeat to advance to the next time step. 



V. RESULTS AND DISCUSSION 

A typical set of simulation parameters are shown in 
Table HI Most of the values are those used by Yamaguchi 



28] 



and Doi 0- According to the Maier-Saupe theory 
of the isotropic-nematic phase transition, the order pa- 
rameter at the first order jump to the nematic phase 
is S = 0.42. The orientational order parameter asso- 
ciated with the polymer backbone is usually a fraction 
of this bare nematic order parameter, which motivates 
our choice of S = 0.3. We should note however that 
much higher values of the order parameter have been re- 
ported for the backbone order parameter in LCEs, par- 
ticularly for main-chain systems where S = 0.9 has been 
observed (2^. The value of r adopted here is lower 
than that used by Yamaguchi et. al. but is consistent 
with creep flow measurements on adhesives manufactured 
from acrylic polymers j30| . 

In Fig. |4ja) we show how the debonding stress cr™ s 
varies with the deformation A (notice the logarithmic 
scale on the x-axis) for three cases; (i) the adhesive is 
isotropic, (ii) the adhesive is nematic with the director 



Atmospheric Pressure (P a tm^ 
Shear Modulus (G) 
Relaxation Time (r) 
Viscosity (r) = Gr) 
Strain Rate (A) 
Surface Tension (7) 
Friction Coefficient (fj.) 
Typical Cavity Radius (Ro) 
Initial Height (H ) 
Length (L ) 
Order Parameter (S) 
Number of Blocks (TV) 



10 5 Pa 
10 s Pa 
30 s 

3 x 10 6 Pa s 
0.1 s" 1 

3 x 10~ 2 J m- 2 

2 x 10 9 Pa s m _1 

10~ 6 m 

10~ 4 m 

5 x 10~ 3 m 

0.3 (0 if isotropic) 

100 



TABLE I. Model Parameters 



initially parallel to the stretch direction and (iii) the ad- 
hesive is nematic with the initial director perpendicular 
to the stretch direction. Each case has the same initial 
random seed for cavity formation. The three cases show 
the typical behaviour of a PSA. We show the behaviour 
up to a very large deformation of A = 20. Whilst such 
large deformations can be achieved experimentally using 
soft, weakly crosslinked polymers [1[ the Gaussian model 
used here does not provide a good description of these 
large deformations, or the detachment of the adhesive. 
However, it does provide a qualitative model of the sub- 
sequent behaviour that may be achieved at lower strains 
in more precise models. We can identify 4 separate be- 
haviours within the plot: 

1. A steep initial rise in the debonding stress until 
A ~ 1.05. 

2. A decrease in the debonding stress until A ~ 2. 

3. An increase in the debonding stress until A ~ 5 for 
the isotropic case and A ~ 10 for nematic cases. 

4. A decrease in the debonding stress for larger 
strains. 

We explain the main processes occurring during each of 
these stages. Fig. 0Ja) shows the local debonding stress 



&lc 



J yy,i 



Pi + P l+ 1 



/Patm, (30) 



the shear strain X xy and the angle 9 that the director 
makes with the y-axis as a function of position in the 
adhesive for A = 1.02. The local debonding stress is 
small near the edges (x ~ ±0.5) of the adhesive and rises 
rapidly towards the centre (x ~ 0). Investigating sepa- 
rately the pressure and polymer stress contributions to 
the local debonding stress reveals the pressure to be the 
predominant cause of the initial rise in the debonding 
stress. As each block expands along the y-direction it 
must also contract along the x-direction. This contrac- 
tion results in large shear strains and concomitantly large 
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-0.5 -0.25 0.25 0.5 -0.3 -0.2 -0.1 0.1 0.2 0.3 



FIG. 4. The debonding stress as a function of the deformation A for the parallel, isotropic and perpendicular arrangements. 
The debonding stress ai oc in units of P a tm, the shear X xy and the angle the director makes with the y-axis as functions of 
position within the adhesive x are shown in (a) for A = 1.02, (b) for A = 1.2 and (c) for A = 3. 



shear stresses. These shear stresses are balanced by the 
pressure difference across the block (see Eq. <JT0j) ) , result- 
ing in negatively large pressures within the bulk of the 
adhesive and consequently large local debonding stresses. 
The angle 9 is largely close to 8 = ±90° for the perpen- 
dicular case - i.e. at A = 1.02 the director has not rotated 
much away from its initial orientation, apart from close 
to the edges x ~ ±0.5. In the parallel case there is ev- 
ident director reorientation, particularly near the edges 
of the block x ~ ±0.5 where the director is at 9 = ^45° 
despite initially being close to 9 = - i.e. towards the 
edges of the block, where the shears are greatest in mag- 
nitude, the director has reoriented to accommodate the 
shear. 



Fig. IDJb) shows how ai oc , X xy and 9 vary with posi- 
tion x for A = 1.2. Comparing with (a) we can see that 
ai oc has become substantially more negative for the per- 
pendicular case while the parallel and isotropic results 
look broadly similar to the situation in (a). The shear 
strains for the isotropic and parallel situation are simi- 
lar in magnitude to (a), but we can now see kinks which 
correspond to cavities within the adhesive. The shear 
strains for the perpendicular case are larger, both when 
compared with (a) and when compared with the paral- 
lcl/isotropic cases in (b). The angle 9 is close to ±90° for 
the perpendicular case while for the parallel case the an- 
gle remains close to 0° around x = but becomes ±50° 
towards the edges x = ±0.5. Fig. [5] shows a representa- 
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FIG. 5. A representation of the blocks for the perpendicu- 
lar (a) and parallel (b) alignments. Some cavities are visible 
within the adhesive layer. The three large circles show close 
ups on the relevant location in the adhesive. The arrows show 
the director orientation. 



tion of the adhesive at A = 1.2 for the perpendicular (a) 
and parallel (b) alignments. Large cavities are evident. 
As can be seen, the shears near the edges of the adhesive 
are larger for the perpendicular case than for the parallel 
case. Note that the traction force on the substrate due 
to the adhesive can become negative under some circum- 
stances. When the shear of the blocks is large, then the 
constitutive model results in the expansion of the block 
and a negative tensile force. The expansion of rubber 
under shear is a well known effect [3l|. It is interesting 
to look at the large central cavity for both the perpen- 
dicular and the parallel alignments. We notice in the 
perpendicular case that the blocks on either side of the 
cavity are sheared in the same sense, while in the parallel 
case the blocks are sheared in opposite senses. This can 
also be seen in the plot of X xy in Fig. QJb), the kink in 
the curves around x ~ 0.05 involves a change in sign for 
the parallel alignment, while there is no change in sign 
for the perpendicular alignment. 

Up until this stage the main contributor to the debond- 
ing stress has been the pressure. As the blocks continue 
to elongate however we anticipate the polymer stress con- 
tribution will become more and more important - for 
a simple nco-Hookean material we would expect these 
terms to scale quadratically with A. The rise in the 
debonding stress after A ~ 2 is principally due to the 
polymeric terms in the stress tensor increasing. Fig.2Jc) 
shows (Tioc, ^xy and 6 as a function of x for A = 3. 
We see now that ai oc is somewhat more uniform as a 
function of position for all three cases. The debonding 
stress is largest for the parallel alignment and lowest for 
the perpendicular arrangement. The debonding stress is 
strictly positive for all positions. It is also interesting to 



note that the arrange of the plot is substantially reduced 
at this strain, indicating that the blocks have slipped. 
For both the parallel and the perpendicular cases we see 
6^0 for all blocks, i.e. the director has largely aligned 
with the stretch direction for all blocks. The shears X xy 
are smaller in magnitude than those presented in (b) , we 
notice however they remain larger in magnitude for the 
perpendicular alignment than for the parallcl/isotropic. 

Above A ~ 5 for the isotropic and A ~ 10 for the ne- 
matic we see a reduction in the debonding stress. This 
is due to the stress relaxation inherent in the constitu- 
tive model of Eq. ([15} . We expect the polymer stress to 
relax over a time-scale set by r — 30s for the isotropic 
situation. In the nematic case two time-scales appear 
t(1 - 5) = 21s and r(l + 25) = 48s. The decay time 
we observe for the perpendicular and parallel alignments 
is set by this larger timescale r(l + 25). At these larger 
strains the director for each block is aligned along the 
y-axis and each block has essentially the same value for 
the local debonding stress (X; oc . 
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FIG. 6. a) The tack energy Sj- as a function of strain for the 
isotropic, parallel and perpendicular cases, b) The relative 
tack of the parallel and perpendicular alignment compared to 
the isotropic case. 

In Fig.[BJa) the scaled tack energy E-y / {AH$P atm ) de- 
fined in Eq. (fTT|) is plotted as a function of deformation 
A for the isotropic, parallel and perpendicular cases. In 
practice the adhesive eventually detaches at a particu- 
lar deformation [27j. As can be seen, the tack energy is 
always largest for the parallel alignment. The tack en- 
ergy is mostly lowest for the perpendicular alignment, 
although this changes when above A ~ 13 where the tack 
energy for the isotropic adhesive becomes lower. This 
plot demonstrates that the adhesive properties of a ne- 
matic adhesive can be altered both by changing the ini- 
tial alignment (i.e. switching from the parallel to the 
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perpendicular geometry) and by changing the order pa- 
rameter (i.e. switching between the isotropic and ne- 
matic phases). The latter transition is relatively simple 
to achieve either by altering the temperature or by us- 
ing photo-active nematics |32j, whilst the former might 
be achievable by mechanically stretching the adhesive 
layer. The difference in tack energy is not so large for 
small deformations, however it can become appreciable 
for larger deformations. In Fig. E]Jb) we show the tack 
energy ratio between the parallel: isotropic and perpen- 
dicular isotropic as a function of the deformation A. For 
the parallel: isotropic plot the ratio is close to unity de- 
formations up to A w 3, after which we see a gradual in- 
crease in the relative tack up to a value of 2.3 by A = 20. 
The curve for the perpendicular:isotropic ratio is slightly 
more complicated. Initially the relative tack is below 
unity up until deformations of around A w 14. beyond 
which the relative tack continues to increase above unity. 
The detachment process is not modelled here so we will 
assume detachment occurs at a strain of A = 10, which is 
consistent with previous work [271 ] . For large strain the 



precise choice of this detachment does not change our 
conclusions. For A = 10 the ratio of the tack energies is 
1.55:1:0.86 for para:iso:perp. 

As can be seen in Table U there are a large number of 
parameters in our model which we might adjust in order 
to maximise the difference in tack between the isotropic 
and nematic states. We now consider how changing 
several of these parameters changes the force extension 
curves and resultant tack. 



A. Varying r 




FIG. 7. Variation of the debonding stress with 
deformation for several values of r shown in the fig 
parameters are as in Table [J 



15 20 



increasing 
ure. Other 



The time constant r is the fundamental time-scale 
over which stress is relaxed away. In the nematic case 
we in fact have two time-scales tj_ = r(l — 5) for the 
relaxation of stress perpendicular to the director and 
Tii = t(1 + 25) for the relaxation of stress parallel to the 
director. Nevertheless the average of these time-scales 
is still t (r = (tii + 2tj_)/3). The time-scale t is also 
related directly to the viscosity r\ which appears in our 
dynamical equation for the cavity radius (see Eq. ([5])) 
via r\ = Gt. Altering r can thus be expected to al- 
ter both the small strain regions of the tack curve where 
the debonding force is largely determined by the cavities, 
and at larger strains where the debonding force is largely 
due to the elastic deformation of fibrils. Fig. [7] shows 
the force extension curves for three different values of r. 
All other parameters are as listed in Table |U We can see 
that the effect of increasing t is to increase the heights 
of the two peaks which occur in the plot and to move the 
peaks to larger deformations. Therefore a larger r tends 
to produce a greater tack energy. Of more interest here 
however perhaps is the relative tack. For A = 10 the tack 
energies are in the ratio 1.57:1:1.11 (||:iso:_L) for t = 10s, 
1.67:1:1 for t = 20s and 1.55:1:0.86 for r = 30s. The 
relative tack values are summarsicd in Table |TTJ 



B. Varying S 

Figure [8] shows the effect of varying the nematic order 
parameter 5 has on the debonding stress as a function of 
the deformation A. All other parameters in these plots 
are as listed in table HI For deformations less than A ~ 2 
there is very little difference the curves for the parallel 
alignment and for the isotropic adhesive. The perpendic- 
ular alignment is quite strongly influenced by the order 
parameter. Increasing the order parameter leads to a 
reduction in the height of the first peak in the debond- 
ing stress. The deformation at which the peak occurs 
also becomes larger with increasing order parameter. At 
larger deformations both the parallel and perpendicular 
alignments are effected by changes in the order parame- 
ter. Increasing the order parameter results in the second 
peak in the debonding curves occurring at a larger defor- 
mation and a larger debonding stress, this can be under- 
stood by considering the time-scale for stress relaxation 
parallel to the director t\\ = r(l + 25). At larger values of 
5 this time-scale is longer, so elastic stresses build up for 
a longer time for increasing 5. For A = 10 the tack en- 
ergies are in the ratio 1.55:1:0.86 for 5 = 0.3, 1.84:1:0.65 
for 5 = 0.5 and 2.07:1:0.41 for 5 = 0.7. Hence there 
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FIG. 8. Variation of the debonding stress with increasing 
deformation for several values of S. Other parameters are as 
in Table U 



FIG. 9. Variation of the debonding stress with increasing 
deformation for several values of G shown on the figure. Other 
parameters are as in Table [J 



is a reversible change in the tack energy by more than 
a factor of 2 for S = 0.7 between the parallel and the 
isotropic states and between the isotropic and perpen- 
dicular states. The relative tack values are summarsicd 
in Table HU 



C. Varying G 

Figure shows the effect that varying the elastic mod- 
ulus G has on the debonding stress as a function of the 
deformation A. All other parameters in these plots are 
as listed in Table |U Note the axes have the same scale on 
all three plots. It is apparent that increasing G increases 
the magnitude of the debonding stress. The deforma- 
tion at which the second peak occurs is the same in all 
three plots, the associated debonding stress at the second 
peak scales proportionally with G. The debonding stress 
at the first peak is larger for larger values of G, however 
the debonding stress at the first peak is not proportional 
to G - the debonding stress at the first peak in fig. [SJc) 
is roughly three times greater than in fig. Ufa) while G 
changes by a factor of 10 between the plots. 

It is apparent that increasing G leads to greater val- 



ues for the absolute tack. The tack values for A = 10 
are in the ratio 1.49:1:0.86 (]|:iso:_L) for G = 3 x 10 4 Pa, 
1.55:1:0.86 for G = 1 x 10 5 Pa and 1.61:1:0.86 for G = 
3 x 10 5 Pa. It is interesting to note that the relative tack 
for switching between the isotropic and the perpendicu- 
lar states is identical for the three values of G listed, this 
is a peculiarity of our choice to quote results for A = 10. 
The relative tack values are summarsied in Table [TTJ 

D. Varying A 

We have previously considered the effect of changing r, 
S and G. These are material parameters of the adhesive 
that can be changed by altering the material chemistry 
or the chain architecture (chain length, entanglement, 
molecular weight, chain branching, etc). We can also 
consider changing experimental parameters such as the 
strain rate A. Figure ITU1 shows the debonding stress as a 
function of the deformation A for several values of A. All 
other parameters are as listed in table |U As can be seen, 
changing A has quite a large effect on the form of the 
curves. For A = O.Ols^ 1 as shown in fig. fTUTa 1 ) there ap- 
pears to be no second peak in the curves for the parallel 
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FIG. 10. Variation of the debonding stress with increasing 
deformation for several values of A shown on the figure. Other 
parameters are as in Tabled] 



and isotropic adhesives, a second peak is just visible in 
the curve for the perpendicular alignment. This form of 
the debonding curve is typical of "liquid-like" debonding 
and is called cohesive debonding. In the current case since 
the Deborah number De = At = 0.3 is less than unity 
the elastic stresses are able to relax away over the time- 
scale associated with the increase in deformation and we 
are left with a viscous liquid type debonding curve. In- 
creasing the strain rate to A = 0.05s -1 produces quite 
a big difference in the curves. The overall magnitude of 
the stresses has increased, and we can now clearly see 
secondary peaks in the debonding stress. The shape of 
this debonding curve is typical of materials intermediate 
between a brittle solid and a viscous liquid, and is known 
as adhesive debonding. Increasing the strain rate further 
to A = 0.1s -1 and we see the general form of the curve 
is the same as in fig. [TOFM, but the overall magnitude of 
the stresses has increased, and the deformation at which 
the second peak occurs is larger - the deformation corre- 
sponding to a relaxation time r is larger for larger strain 
rates. These observations are consistent with experimen- 
tal results on isotropic adhesives [33] . It is clear that we 
obtain larger absolute tack values for larger strain rates. 
Comparing the relative tack values at A = 10 we ob- 
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0.86 
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0.86 
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1.55 


0.86 
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1.61 


0.86 




0.3 


1.49 


0.86 


A( s -1 ) 


1 


1.55 


0.86 




3 


1.61 


0.86 



TABLE II. The relative tack values at a deformation of A = 10 
for parallekisotropic and perpendiculanisotropic for various 
parameters. 



tain 1.23:1:1.03 for A = 0.01s -1 (||:iso:_L), 1.72:1:1.09 for 
A = 0.05s -1 and 1.55:1:0.86 for A = 0.1s -1 . The relative 
tack values are summarsied in Table [TT] 



VI. CONCLUSIONS 

In this paper we have employed a modified version 
of the block model of Yamaguchi et al. Q and the ne- 
matic dumbbell constitutive equation of Maffctone and 
Marucci [2l[ to model the adhesive debonding of nematic 
elastomer adhesives. These models include several ap- 
proximations and simplifying assumptions in their de- 
scription of adhesive debonding: 

1. The flow of films is assumed to be a superposition 
of slippage and a piecewise linear deformation of 
the block, 

2. The cavities are assumed to be described by gaps 
between adjacent blocks, 

3. The cavity dynamics are modelled by the Raylcigh- 
Plesset equation, 

4. The slip velocity at the substrates is assumed to be 
a linear function of the shear stress, 

5. The nematic director reorients very quickly com- 
pared to the polymer relaxation time, 

6. The polymers can be described as Gaussian chains, 

7. The adhesive cannot debond from the surface or 
rupture. 

Despite these simplifications the main physical processes 
which occur during debonding are captured and the over- 
all result of our modelling, i.e. demonstrating a reversible 
difference between the tack energy as a result of the ne- 
matic to isotropic transition is valid qualitatively. 
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Within these modelling assumptions we have shown 
that there can be substantial differences in the tack en- 
ergy when comparing an isotropic and a nematic adhesive 
which are due to the difference in their rheology. For ex- 
ample with a nematic order parameter of S = 0.7 we were 
able to achieve a relative tack energy at A = 10 of 2.07:1 
between the parallel aligned nematic and the isotropic 
and 0.41:1 between the perpendicular aligned nematic 
and the isotropic (and thus a ratio of 5.05:1 between the 
parallel and perpendicular alignments). Given the ease 
with which one can reversibly cycle between an isotropic 
and nematic phase this work gives impetus to the ex- 



perimental investigation of this mechanism of switching 
adhesives on and off. 
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